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Abstract 

We extend the low-energy effective field theory relevant for & — > s transitions up to operators of mass-dimension 
7 and compute the associated anomalous-dimension matrix. We then compare our findings to the known results 
for dimension 6 operators and derive a solution for the renormalization group equations involving operators of 
dimension 7. We finally apply our analysis to a particularly simple case where the Standard Model is extended 
by an electroweak-magnetic operator and consider limits on this scenario from the decays Bs — > l-i-^ fJ-' and 
B Kuv. 
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1 Introduction 



Flavour-violating processes are well-known as a central test of the Standard Model (SM). The pattern conceded 
to flavour transitions is indeed particularly constrained in this model, the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix encoding the only source of flavour-breaking effects while only the charged currents of the weak interaction 
convey flavour-violation at tree-level. Consequently, flavour transitions also deflne closely watched observables 
in the quest for new physics and, due to the absence of positive deviations from the SM predictions, set serious 
constraints on the forms that physics beyond the SM (BSM) could take. One of the latest results is the observation 
by the LHCb collaboration of the decay -> fJ-^fJ-^ |1 , with a branching ratio very compatible with the SM 
expectations (refer to [2j for a recent summary): 

BR{B° ^ = (3.23 ± 0.27) • 10"^ ; BR{B° ^ M+Z^")^^ = (3-56 ± 0.18) • lO'^ (1) 

On the theoretical side, dedicated tools have been devised in order to study flavour- violating processes, in the form 
of low-energy effective field theories (EFT; refer to e.g. [3 for a review). These EFT's allow for a separation of the 
long distance, low-energy strong interaction effects and the short-distance, flavour-changing physics. In the context 
of S-physics, at the level of terms of dimension 4 and smaller, the relevant EFT consists in a QEDxQCD model 
with flve quark-flavours (u, d, c, s, b) and 3 charged leptons (e, /i, r), as well as neutrinos. Then, in the low-energy 
processes involving those flelds, at a scale fib ~ Mb, the impact of higher-energy (top/Electroweak/Higgs/BSM; 
we assume here that there are no further low-energy 'invisible' flelds) physics can be essentially encoded within 
operators of dimension > 4, collectively deflning the 'effective hamiltonian'. The strong-interaction problem then 
consists in evaluating the 5- matrix elements driven by these operators among physical (mesonic/baryonic) states: 
this question is answered, either through lattice-QCD, QCD sum rules, heavy-quark expansions and other theoreti- 
cal descriptions of the non-perturbative strong-interaction effects, or phenomenologically, through an identiflcation 
of the decay-constants by comparison with a few standard channels. The short-distance problem is summarized 
within the couplings multiplying the operators. Those must be matched at high-energy with the predictions of the 
'more fundamental' theory (Standard Model, supersymmetry-inspired models, etc.): scattering amplitudes in both 
the EFT and the 'full-theory' are equated at the matching scale ^ M\y, hence deflning a boundary condition 
for the parameters of the EFT in terms of those of the underlying model. 

To relate these two scales, /ib and /ioi one relies on the Renormalization Group Equations (RGE) driven by the 
renormalization of the EFT: leading logarithmic contributions can thus be consistently resummed. It is mostly in 
this part of the procedure that we will be interested in the following. 
To flx notations, let us write the lagrangian density of the EFT under consideration: 

Ceff = -If.^F'^'^ - ^G^.G'^^^ + ^fb''D, -inf]f ^ nf^f^ (2) 

with f — u,d, c, s, 6, e, /j, r, Ve.fi.r, — d^ — leQfA^^ — «.gsGJ^ the covariant derivative, ruf the mass of the fermion 
/, Qf, its charge; F^, = d^A, - d,A^ and G^, = G^^T- = (9^GJ - d^G^ + gsr^'G^GDr- represent the 
electromagnetic and gluonic fleld tensors, respectively; T" and /"'"^ denote the SU{3)c generators and structure 

2 

constants; e and gs respectively stand for the elementary electric charge and the strong coupling constant; a = 
as = ^] "Hgj™"^^ is the effective hamiltonian. 

The question arising at this point is that of the operators that one needs to consider in order to describe 6 — s 
transitions. Sensibly, people have considered, up to now, only operators of the lowest possible mass-dimension, 
that is dimension-6 operatorf|^ This approach is justified by the suppression factor of ^ 5 • 10~^ which is 
expected for higher-dimensional operators. Moreover, it was (with reason) regarded as sufficient to confine to the 
smallest subset of dimension-6 operators that would close under renormalization and for which the classical models 
(SM, supersymmetry-inspired, etc.) would generate a non-trivial contribution: 

<7^'^-E^^(/^)o1"^(a') (3) 

i 

where the hst of operators Oi can be read in e.g. [s], [i] or [s] (with small variations; note that a factor Gf/V2 
is conventionally factored out in the usual notations). The determination at leading order of the anomalous- 
dimension matrix for the four-quark operators of this subset is quite old: [6 - 10 . The renormalization of the 



magnetic and chromomagnetic operators can be found in |11 , while 12 included the mixing of those with the 



^In fact, the magnetic and chromo-magnetic operators have mass-dimension 5. However, they always appear, e.g. in the SM, with 
an additional mass-suppression, lowering their scale to an apparent dimension 6. 
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four-quark operators (a two- loop effect). One may refer to [13 ■ 15 for the analysis of the semi-leptonic operators. 
Later works have focussed on next-to-leading order (0(as)) [5 16-20 , electroweak (0(a); see e.g. 21 ) and finally 



next-to-next-to-leading QCD (0(a|); refer e.g. to the summary in 22 ) effects: this formidable amout of work 
allows for a theoretical prediction, e.g. in the SM B — ;> X^j decay, competitive with experimental bounds. 
In this paper we choose to adopt a different, more unprejudiced if somewhat more anecdotical, approach to the 
renormalization of the 6 — > s EFT: we shall consider all possible (on-shell) operators up to mass-dimension 7 
and compute the corresponding anomalous-dimension matrix at one-loop QCD order. As far as we know, no 
attempt has ever been made in that direction, due to the suppression of the order ^ 5 • 10~^ which one 
expects for higher-dimensional operators. In fact, if one considers the matching conditions in the particular case 
of the SM, with its restricted flavour-changing currents, the suppression would be even larger. The inclusion of 
higher-dimension operators hence admitedly appears in this concrete case more as a curiosity than a compelling 
necessity. There are however several reasons why analysing dimension 7 effects may not be completely irrelevant: 



1. From the point of view of precision physics, one observes that as{Mz) 



jj^: with increasing precision in 
the SM evaluation as well as experimental measurements, observables shall eventually become sensitive to 
dimension 7 effects. 



Certain new-physics contributions are actually 'hidden' dimension 7 effects; a simple example lies in the fa- 
mous Higgs-penguin contributions to dimension 6 &sZZ-operators 23 , relevant e.g. in supersymmetric models 
at large tan/3: the corresponding coefficients actually contain a factor to;,, formally increasing their order 
to dimension 7. Note however that other dimension 7 effects are not expected to receive an equally large 



tan/3-enhancement, although they should be relevant already at subleading order 24 



New-physics in dimension 6 operators is already stringently constrained, essentially enforcing the usual 
'Minimal Flavour Violation' condition. A possible strategy to account for this absence of new-physics in 
flavour observables would be to reject it on operators of higher dimension. Similar proposals have been 
made in the neutrino sector to concile neutrino masses, baryon/lepton- number and Icpton-fiavour violating 
Note that, here, we do not propose a mechanism that would ensure the suppression of new 
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processes 

physics by rejecting it on operators of dimension > 7 (although this might be achievable, e.g. by assigning 
adequate charges under a discrete symmetry), but we simply mention this possibility as a motivation to 
consider such operators. 



4. If one parametrizes new physics blindly in an expansion of SM dimension 6 operators 26|27 , it turns out that 
certain operators would only have a dimension-7 signature at low energy. Including dimension 7 operators 
for the 6 — >■ s transition thus naturally enters an unprejudiced analysis of physics BSM. 

In the next section, we shall derive a list of all the (on-shell) operators of mass-dimension 5, 6 and 7 which may 
intervene in the 5 — > s EFT. We shall then detail the calculation of their ultraviolet (UV)-divergences at one-loop 
QCD order, before we proceed to the renormalization and establish the RGE's. The following section will be 
dedicated to the solution of these RGE's in some specific cases. Finally, we shall illustrate our discussion by 
presenting a concrete, if naive, case where our analysis of dimension 7 operators apply. A short conclusion will 
eventually summarize our achievements. 



2 Operators of dimension 5 — 7 in the 6 — > s transition 
2.1 List of operators 

The first step of our analysis consists in establishing a list of all the operators of dimension 5, 6 and 7 intervening 
in the 5 — s transition (but Lorentz + gauge invariant!). For simplicity, all the fields shall be taken on-shell, 
i.e. satisfy their equation of motion: this will be sufficient, at least for the leading-order calculation that we aim 
at. For a discussion concerning the relevance of on-shell EFT's and in particular the question of applying only 



'naive' classical equations of motions (dismissing ghosts and gauge-fixing terms), we refer the reader to 28 . We 
can obviously distinguish among three categories of operators: four-quark, two-quark -I- two lepton (semi-leptonic) 
and 6s-Gauge operators. 

Let us start with four-fermion operators: the fermion fields already account for a mass-dimension 6, leaving room 
for at most one covariant (for gauge-invariance) derivative, when one restricts to operators of mass-dimension < 7. 
Considering chiral fermions (that is, we include projectors Pl,r = ^"^^ in the fermion products), there are at 
most three ways to contract the spinor algebra: 1 provides scalar currents, 7'^, vector currents and a'^'^ = [7^,7''], 
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tensor current^ any higher combination of 7-matrices can be reduced down to those three cases (through the use 
of the Levi-Civita tensor e^i/po- and identities of the Dirac algebra): note that 75 only gives a sign, when applied 
on chiral fcrmions. Moreover, we may always keep the b and s of the flavour transition within the same current: 
other combinations are made redundant by the Fierz identities (see e.g. Appendix [A|) . 
Those considerations allow us to construct the relevant three classes of dimension 6 operators: 

1. scalar {hPL^Rs){JPLMf)\ 

2. vector {h^' PL,Rs){h^,PL,BJY 

3. tensor {ba^'^ PL,Rs){fafj.,^PL,Rf)- Note that only two chiral combinations are possible for the tensor currents: 
L X L and R x R; the other combinations, L x R and R x L, are identically zero. 

One may then consider authentic dimension 7 operators by incorporating one covariant derivative in the fermion 
products. Two possibilities appear: 

1. contracting the Lorentz index of with a vector current 

2. contracting it with a tensor current, the second tensorial index being contracted with the second, vector 
current. 

Note however that, using the equations of motion on fermions {{tP — rrif)/ = 0) and realizing partial integrations 
(i.e. adding a total derivative to the lagrangian density), the resulting set of operators is largely redundant (together 
with the dimension 6 operators). Indeed, the second possibility which we mentioned {'D^j'^ ^ct^lu ) can always be 
reduced down to operators of the first class ('D^ ® 7^^') plus dimension 6 operators (multiplying fermion masses). 
Additionally, certain combinations of the operators of the first class reduce to dimension 6 terms. We thus retain 
only two kinds of linearly independant operators: 

1. [bl(d - 'DYPLM{h^.PL,Rl)\ 

2. {h''PL,Rs)[HD - 'P)^,PL,Rf]■ 

For semi-leptonic operators, i.e. when / is a lepton (f = I), the analysis is essentially over. 

Let us therefore focus on four-quark operators {f — q). One should then also consider the contraction of the colour 
indices. Two possibilities arise: 

1. product of two colour-singlet currents; colour-indices contracted between the b and s on one side, q and q 
on the other; 

2. product of two octet currents: colour-indices contracted between the b and q on one side, q and s on the 
other. 

In the special cases where q = b, s, however, Fierz identities make this distinction superfluous, so that we may 
consider only singlet products then (refer to Appendix [A|. This is our final word for four-fermion operators. 
We now turn to &s-Gauge operators. The methodology follows that of |12j for the dimension 6 operators: one 
may simply write all the possibilities to include covariant derivatives (at most four for operators of dimension 
< 7) within the fermionic current. Using partial integration and equations of motion, it turns out that all these 
covariant derivatives can be combined in field-strength tensors. One thus simply needs to consider the possibilities 
to combine the indices of these field-strength tensors with those of the fermionic current: 

• When only one field-strength is present, it can only contract with a (Lorentz) tensor current, and, in the 
case of the QCD field strength C^^, with a SU{3)c octet current. 

• When two field-strengths are present, we may either contract their Lorentz indices together - thus reducing 
the fermionic current to a scalar - via the metric or a e^i/po- tensor -, or contract two of their Lorentz indices 
with a tensor current (the other two through the metric, or equivalently e^vpa)- In the case where only 
one QCD field-strength is involved (among the two), one needs again a SU{2i)c octet on the fermionic side. 
When two QCD field-strengths are involved (G^^Gpg.), the color indices may be contracted together {5°''^: 
colour-singlet), with a symmetric tensor ({T°, T''}) or with an antisymmetric tensor ([T°, T'']) (superposition 
of colour-octets -I- singlet), depending of the compatibility of these structures with the Lorentz form of the 
fermionic current. 

*Note that the definition of cr^"' which we adopt here for simpUcity differs somewhat, by a factor r/2, from the one the reader may 
have encountered in the Uterature. 
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We now have all the ingredients to present the list of relevant operators. Note that their normalization is a priori 
free: our choice can be justified a posteriori to ensure that all the RGE's intervene at the same, leading-order in 
as- Note however that this choice may be slightly misleading, for instance in the SM, for reasons that we will 
discuss in the next section, when we solve the RGE's. 



{bs){ll) operator^ 



SlRL,R^fzrn,{bPLMmPL.Rl) 



^L.R L,R - 
^L.R L,R - 
^LM L.R 
fjl 

^LM L.R 



f^mtibrPL.,Rs){h^.PLMl) 

f^mb(ba>'''PL,Rs){la^,PL^Rl) 



(4) 



:{h''PLMs) 



I 



il5+2QieA] Pl.rI 



{bs){qq) operators 



C9 

'-'l.r l,r 

^L,R L 



= ■mb{bPL,Rs){qPL.Rq) 
^ = m.b(b-i^'PLMs){qiy,PLM.q) 
'Plr LR^ '"T'biba'"' PL,Rs){qaf,^PLMq) 



^LM L,R 
^LM L,R 



Pl,rs 



ih^PL.Rs) 



iqip.PL,Rq) 



qi 



PL,Rq 



^L,R L,R = '^b(haPL,RSl3){qfiPL,Rqa) 
^L,R L,R = 'm'bibaj''PL,RSl3){qpJ^PL,Rqa) 

= mb{baaf^'^PLM.sp){qi3cr^^PL,Rqa) 



'L,R L,R 
1/1 

'^L,R. LM 
njq 

^L,R L,R 



D] P, 



{qpi^iPh.Rqa) 



[Ki^Pl.,rsp) 



PL.Rqo 



(5) 



(6s)-Gauge operators: 



E 



L,R 



Ql,r - ^ 

cS _ _a " 
'^L,R. 



^L,R 

F'^ - 
^L,R - 

rf-L.R 

njS 
"■L,R 

T~l-L,R 



is 



L,R 



"""" ha'"' Pt j}sF 
™' Io^^'^T'^Pl^rsGI^ 



S _ , a 



f-^hPL,RsF^,F^^- 
^bPL,RsF^,F>^'^ 
f^ba-pPL^RsF,,F^^ 
-- ^bT^PL^RsF^.G^t^'^ 
- fLlT-PL^RsF^.G'^f''' 
. f-ba-^PT-PL.RsF'iG'^.p 

QL,R^bPL,RsGl,G^>^^ 

Ql,R^^bPL.,RsGl,G-^^'^ 
QIr = bi^^PL,RsG?..,G''^'' 



QlR^^b^^^PL.RsGl,,G'>^^- 



(6) 



This list determines the effective hamiltonian of our EFT (the generic notation O* spans the whole list): 



n 



(di.„<7) 
eff 



^a(/i)0^(A*) + /i.c. 



(7) 



^We repeat that, for the tensor operators, only the L X L and R X R combinations are relevant, which may not be obvious from 
our notation. 



^We discard S,V,T,'H,'H for q 



, s since they reduce to S, V, T, H, H'^ due to Fierz transformations which are provided in 



Appendix K| The colour indices o, /3 of the fermions are displayed when not trivially contracted. Note finally that the ambiguous 
notation ^ D y,qa stands actually for {D^q)a. ■ ■ 
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Summing over all the flavours and chiralities, we count 251 terms. 

For the purpose of illustration, let us sketch the typical matching conditions that one could expect for these 
operators in a Minimal-Flavour- Violating model (the CKM matrix elements are written as Vqq'). Only the op- 
erators (bs){qq) with q — c,u would receive a contribution at tree-level: Cq=c,u ^ '^^m^^'''' (^^^ '^^^^ q = u 
may be neglected because of the CKM suppression). All the other contributions would arise at the loop- level, 
Cdim 6 ^ '^^m'b^"' ^'^'^ ^dim 7 ^ for thc opcrators present at dimension 6 and 7 respectively. 

2.2 UV-divergent amplitudes involving the dimension 5 — 7 operators 

Before renormalizing the EFT, we compute the divergences associated with QCD loops. We perform this calcula- 
tion in the most naive conceivable way: in the Feynman-t'Hooft gauge and dimensional regularization D = 4 — 2e; 
we keep only the divergent terms [2 — y] . Note that we will assume that the couplings Ci come together with 
an electric charge factor e whenever a photon or a lepton appear in the external lines: in this fashion [bs){ll) shall 
be regarded as electroweakly suppressed with respect to (bs){qq), rejecting (bs){ll) UV-contributions to (bs){qq) 
operators to a subdominant order. The computation can be organized in 'blocks': 

• (bs){ll) contributions to the {bs)(Jl) amplitude: see Figjl] 




Figure 1: Contributions from the (bs){ll) operators to the bs — )■ U amplitude. Here as in the following plots, the 
grey blob represents a vertex associated with the operators of mass-dimension 5 — 7. The subscript 'H\ i/'' under 
some of the diagrams indicates that such diagrams are relevant only for this kind of operators. Feynman diagrams 
are plotted with JaxoDraw . 



• (bs){qq) contributions to the (6s)-Gauge amplitude. 




a 



Figure 2: Contributions from the {bs){qq) operators to the bs — 7/(7 amplitude. Again, as in the following figures, 
the label ' iJ*^, g'?, H*^, H'^' indicates that the corresponding diagrams intervene for such operators only. 

Such contributions involve a quark loop. While the external fcrmions are explicitly taken on-shell, through 
the application of their equations of motion, explicit off-shell terms appear for the external photons and 
gluons: such terms must organize as an equation of motion for the gauge-boson line, which provides us with 
a cross-check of our calculation. They do not matter for the renormalization so that we will not keep them 
explicitly in the divergent amplitudes. The amplitudes (bs) — 7* or (bs) — g* are relevant for the divergent 
amplitudes {bs){ll), (bs){q'q') and (bs){gg) though, which we will present later. 
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Figure 3: Contributions from the {bs){qq) operators to the bs 2j amplitude. +{fi O v) signals that one must 
add the 'mirror' diagrams exchanging the two photon lines. 




Figure 4: Contributions from the {bs){qq) operators to the 6s — »■ 7ff amplitude. 
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+{fi, a ^ u, b) 



9 a 





Figure 5: Contributions from the {bs){qq) operators to the bs — > 2g amplitude. 



The contributions to the amphtudes bs 7/5 are shown in Figj2j The diagrams intervening in bs 
27/73/23 are depicted in Fig |3|4|5 Note that in the case of the dimension 6 four-quark operators, no 
contribution to the bs — VV (where V, V G {7, g}) operators is expected: the dimension 6 basis is stable 
and does not require dimension 7 counterterms. The corresponding amplitudes in Fig |3|4|5| therefore provide 
a simple check of the matching conditions determined in Fig[2j For the authentic dimension 7 operators, 
however, the calculation of the bs — VV amplitudes is fully relevant. 

The operators {bs){qq) with q = b, s deserve a particular attention. Beyond the 's-channel' diagrams, similar 
to those of the other (bs){qq) operators, they indeed allow for a 't-channel' contribution: these can be 
viewed as '^-rotated' contributions associated with the S'^, V, T^, "W, H'^ operators, which, in the case of 
q = b,s, coincide with S'', , T'', H'', H'^, via Fierz identities. Instead of computing the new 't-channel' 
contributions]^ we thus use our generic result for the 5', V, T', 'W, H'' operators and project it on S"^, , 
T'^, H'^ , H'^ with the Fierz matrices ^pjgj.^ exchanging the two bases and defined in Appendix |Aj A factor 
— 1 appears in this process corresponding to the anticommutation of two fermion fields / the transformation 
of an open fermion loop into a closed one. 

(bs)(qq) contributions to the {bs){ff) amplitude. 

The contributions to the bs — >■ II amplitude proceed only from the exchange of an off-shell photon, as 
depicted on the first diagram of Figj6j For the bs q'q' amplitude, contributions originate similarly from 
the exchange of an off-shell gluon (note that the photon exchange would be of a higher-order). Additional 
'diagonal' contributions are obtained for q = q' through the dressing of the {bs){qq) vertex by a gluon line. 
All the corresponding diagrams are depicted in Figj6] 



''' We, in fact, checked this relation explicitly. 
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.5- 




q s. 



Figure 6: Contributions from the (bs){qq) operators to the bs — )■ II and bs — > qq amphtudes. The grey blob with 
a quark loop stands for all possible (bs){qq) contributions to an off-sheU photon/gluon: refer to Figjg] 







Figure 7: Contributions from the (6s)-Gauge operators {E, Q and H O 71^ to the bs ^ ^ amplitude. We 
also define the b ^ s self-energy associated with the operators Q (and Q ^ , QP ■, , , Q^) and depicted 
on this figure, as well as the following ones, by a fermion line with a grey blob and a gluon bubble. 
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• (6s)-Gaugc contributions to the (6.s)-Gauge amplitudes. 

Once again, off-shell photons and gluons appear within the calculation of these amplitudes, offering a nice 
crosscheck of the result, given that they must combine into an equation of motion for the corresponding 
field. Although such terms will not matter for the rcnormalization, bs — > 7* and bs — > g* will be relevant 
for bs — >■ ll/qq/gg- Another check, particularly in the case where two gauge bosons appear in the final state, 
is simply the projectibility of the result on the basis of operators that we have defined: indeed the results 
typically reach this form only after the summation of several diagrams, the coefficients and colour-factors of 
which must combine in the appropriate way. Note that only 'authentic' dimention 7 operators need to be 
included in the amplitudes bs — )■ 77/75/95: the presence oi E ov Q in the divergences of such amplitudes 
(substracting however the countcrtcrms of E and Q) would signal an instability of the dimension 6 basis, 
which would require dimension 7 counterterms. This evidently does not occur. 




Figure 8: Contributions from the (65)-Gauge operators {Q and Q) to the bs ^ g amplitude. 




<-»■ v) 



Figure 9: Contributions from the (6s)-Gauge operators {£ ^ S^,£^,£'^ and H) to the bs 27 amplitude. The 
black blobs signal a counterterm from the bs^ or bsg amplitudes. 
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Figure 10: Contributions from the (6s)-Gauge operators {H and Q) to the bs 75 amphtude. The (. . .) stand 
for additional diagrams involving a gluon tadpole, which give vanishing contributions. 
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Figure 11: Contributions from the (5s)-Gauge operators (Q) to the 6s — > 2g amphtude. The grey blob with an 
inscribed gluon bubble and attached to a photon/gluon line stands for all the contributions to ofF-shell photon/gluon 
production: refer to FigfZlSl 
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• (5s)-Gauge contributions to the (6s)(//) amplitudes. Such contributions proceed from off-shell photon or 
gluon exchanges or a gluon 'pseudo-boxj^ and are depicted on Fig 12 




Q,Q 





Figure 12: Contributions from the (6s)-Gauge operators (Q, H and Q) to the bs — » II /qq amplitude. 

Finally, we may cast all these divergent (on-shell) amplitudes into the following matrix form: 

where (C) = {Cgl^^'^ , . . .) is a row vectoJ^ collecting all the couplings of the effective Hamiltonian and (F) is a 
column vector constituted by all the form-factors associated to the operators and intervening in the amplitudes 
that we have just presented (we only need one form-factor per operator since we have not considered redundant 
amplitudes such as bs -> //7, etc.). M depends only on fcrmion masses, gauge couplings and charges. 

2.3 Renormalization of the operators 

2.3.1 Basics of the renormalization of the EFT in the MS-scheme 

Let us recall a few basic ingredient of the renormalization procedure for our EFT in the MS-scheme. 
Dismissing the effective hamiltonian for the time being, the (dim< 4) EFT is a renormalizable QEDxQCD model. 
We neglect the QED loop-effects here, since a <^ as, hence focus on a QCD renormalization only. The fermion 
wave-function and mass renormalization constants, Zf and Zmf respectively, are determined by requiring that the 
counterterm contribution to the / self-energy (with p external momentum), Y]'^"^{p) = i(Z f — \)i> — i{Z fZ^j — 1), 
cancels the divergent loop contributions: I]'r'"'"'(p) = (leptons) or _ »C'2(3)as _|_ 4^^j (quarks), where 

C2(3) — I denotes the Casimir operator of SU{?>)c in the fundamental representation: 

The gluon self energy (involving gluons, quarks and ghosts) similarly determines the gluon renormalization con- 
stant: in the Feynman gauge, ~ 1 ^ ^ \%np — |-/Vc], with np = 5 the number of quark flavours and 

Nc = 3 the number of colours. From the quark-gluon vertex, one renormalizes the strong coupling constant gs, 

2 

resulting in the RGE (from the scale-independence of the bare coupling) for as = with /x the renormalization 
scale: 

das nfn D\ . o, \ Of \ 2/3o 2 , ^/ 3^ o llNc-2nF 23 



2 2-- a5+/3(as) ; p{as) = -^a's + 0{a's) ; /3o 



«^(^)^ (10) 



*By which we mean the last diagram of Fig|l2l 

®We refrain from writing this row vector witti a transposition symbol (i.e. (C)^), not to make already-heavy notations as will 
appear in the following sections completely unreadable. 

^"Those are well-known basics of the QCD renormalization and are available in countless textbooks. We mention them only for the 
sake of completeness as well as clarity in our notations. 
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The bare fields Zff as well as the bare masses Z^fiTif must be scale independent (note that the /^-dependence 



enters Zf and Zm^ indirectly through as)- One deduces the quark-mass running: 



dmq as 



0(4) ^ 



drUq 
din as 



3C2(3) 

/3o 



mq + 0{as) 



3C2(3) 



(11) 



Let us come back to the effective hamiltonian. The loop-divergences in the amplitudes, Eq.Q, must be compen- 
sated by the operator counterterm contributions: i{C)[Z'~^ — l]Zr{T), with Z"-^ the operator-renormalization matrix 
and (the diagonal) Z^ the renormalization constant of the fields (and normalization) entering the form factor T. 



c _ 



[1 



M 



This determines: Z' ~ \m. 7 — „ 

4,r(2-f 

ZpOr is scale-independent) leads to the RGE for the couplings (C): 



]Zp ^. Scale-independence of the bare couplings {C)Z^ (by definition of Zr, 



d{C) 
din /I 



{C)ic 



d\Z^ 



ic 



din II 



C\-l 



= [1 - 
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47r 2- 



(12) 



where the /(-dependence of Z'-^ originates from the dependence of A4 on gs and quark-masses, as well as from the 
renormalization constants in Zr- We have introduced the anomalous dimension matrix 7c. 

Considering the normalization of the operators that we introduced at the begining of this section, one may factor 
out harmoniously C2{S)as in the matrix M, leading to: 



div — 



iC2i3)as 
f) 



An 2 



{C)M{r) 



(13) 



where is a constant matrix (depending only on quark-mass ratios ^)- 

At this leading order in QCD, using EqjlOj Eq|T2] may be written explicitly as: 



d{C) ^ 2agC2(3) 
d In /i 47r 



(C) 



V-M 



0{al) ^ 



ic 



djC) ^ C2(3) 
din as /3o 
2asC2(3) 



(C) 



where we have introduced the diagonal matrix ^"sC2(3) -p 



V-M 



dZy 



V-M 
+ 0(a|) 



0{as) 



(14) 



^ . The corresponding scaling may be extracted 

from the field renormalization factors (Zq^'^ for each external quark and ^3^^ for the external gluons) and the 
normalization factors (inserting gs, as, Tnt). 



2.3.2 Renormalization Group Equations for the operator basis 

We can now explicitly extract the RGE's for the operators under consideration. The operators are ordered as 
above in Eq.( 4|5|6 ). Moreover, as far as the chiralities are concerned, we use the ordering LL, LR, RR, RL for 
four-fermion operators {LL, RR for the tensors) and L, R for the 6s-Gauge operators: we will not write them down 
explicitly in the following, so as to keep notations as tractable as possible, but remember that the coefficients 
come together with one or two chirality indices. The anomalous dimension matrix can be split in several blocks, 
corresponding to the various types (four-quark, semi-leptonic, 6s-gauge) of operators: 
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2asC2(3) 
An 



V" -M" ■ 

-^M9l -M^i V99 - M^9 



(C) = (C',C"?,C3) 



(15) 



where the meaning of the subindices I, q and g should be transparent. For the sake of clarity, we choose to present 
the various blocks along with the specific RGE that they affect, separating the diagonal scaling contributions from 



the UV-divergent diagrams of section 2.2 However, should one wish to implement the whole 251 x 251 7c matrix 
directly without considering our detailed study, the corresponding results are gathered in Appendix [B) 



{bs){U) couplingj"^ 



^Noto that operators with different lepton flavours do not mix, justifying the use of only one index I. 
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The matrices M'-^^siM must be extracted from om calculation of the divergent contributions^ 
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We define the object [1; — Vpjgj.^](q=h_s) as follows: when q 7^ 6, s it is simply the (36 x 36) identity in the space 
defined by the corresponding {bs){qq) operators; when q — b,s, then it is a (18 x 36) matrix built with the (18 x 18) 
identity in the left-hand subblock and the (18 x 18) transition-matrix — Vpj defined in Appendix AI in the right- 
hand subblock. Its origin is related to the facts that, for q = b, s, Fierz identities shorten the list of independent 
operators (compared to the cases q 7^ 6, s), and that the additional diagrams 'in the t-channel', which open for 
q = b, s, can be related to contributions that the redundant (omited) operators would have 'in the s-channel': refer 
to our discussion in section [2.21 
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The objects [l,0](g=f„s) and 



are again defined as the (36 x 36) identity when q ^ b,s and the (18 x 36 

iq=b.s) 

and 36 x 18 respectively) matrices defined by the two (18 x 18) subblocks corresponding to the identity and the 



^■^To make our notations more transparent: each entry in the matrices JVi corresponds to a4x 4, 4x2, 2x4 or 2x2 block in 
chiraHty space corresponding to the type of coefficient it multiplies (row index) and the one it affects (column index). These blocks 
involve the 4 X 4 or 2 X 2 identities, 1, as well as several other matrices, S. 5 . <I>, etc., which we define alongside their first site of 
appearence. When we mention operators alongside the matrix Ai^' , in Eq. jl8| , it is simply in order to clarify which type of operator 
corresponds to the given row (since we cut the matrix to skip numerous entries). 

^^Note that operators with different quark fiavours mix, justifying the use of two quark indices q and q' . 
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null matrix, when q = b,s: their meaning should be clear considering the shortened base in the case q = b,s, due 
to Fierz identities. Prom our calculation of the divergent contributions, the matrices A^' read as: 



q—b,s 



+ [l-^Fierz]<^'=M^^lt 



^ierz 



(21) 





81 




u 


1 A 


u 


u 


U 







3 A 


u 


u 









O-Tl 1 3 \ ^ 

21 + 4^-3 


U 


u 


u 


U 




u 


u 


u 




24A^ 














-72A^ 


















Ah 




n 
U 


n 
1 


U 




Dh 




U 


u 


u 


Mi!, = 


Ah 

fl 


R - 
^h 

n 
u 







u 

U 


H 

u 


5-11 
4I 


Dh 








n 

n 
u 


n 

u 









8 V-^ ^ zj3; 














^(431-211:3) 













-36A^ 





-fl 








-84A^ 







fl 










An 


B-H 













Dh 







1 







. An 

























1 




"0 

























0" 







































































































univ 


1 

































~ 4 






































1(1 + s) 













-(1 + E) 






















































1(1 + S) 








- 


(1 + S) 







































Ah^ 



Ah 



A-H 



B 



H 



Bu 



1 

4 



An = ~ 



3 



3(1 + E) + S3 + 
31 + S3 + ^(3S + ^) 

3(l + S) + S3+£; 

3^(3S+"£) 
81 + S3 + ^(-8S + ^) 
-8(1 + S) + S3 + 



3 

8 mb 

1(31 + S3) 



(31 + S3)-|^(3S + i;) 



Ch = - 



' - = 

Dh = 
Dh^ 
Dh = 

Dn = 



3 

8 m;, 
3 

4 mfc 

3 



7 m, 

8 mb 

-1(31 



3(1 + S) + S3 + £ 
3(1 + S) + S3 + 

(31 + S3) + ^(3S + i;) 
3(1 + S) + S3 + E 
f S3 



E) 



|^(3S + ^) 



-!(S3 4 

1(111 -7S3) + ^(llS-7i;) 
^ [ll(l + S)-7(S3 + E) 



S3 = 



E = 



10 

0-100 
10 




A = 



1 




1 





1 
-1 




1 



1 





1 







-1 





-1 







15 



rrib 

nib 
11 rrig ^ 

8 nib ~^ 
11 niq 6 
4 mfc " 









12 V" ^ ^7 
1 _|_ 



6 ■ inb ' 












64 nib 



U 


U 


U 



































15 "tg 


8 nib 








15 m,^ 


4 nib ^ 












1-10 
0-11 



The definition and meaning of the notation 
• (6s)-Gauge couplings: 
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should be obvious by now. 
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The matrices ^ss-^f proceed from our calculation of the divergent contributions: 
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= [1; 



•^FierzJ 



























































































(25) 



3 Solving the RGE's 

Having presented, in the previous section, what is meant to be the main result of this paper, i.e. the RGE's for 
all operators of mass-dimension < 7 intervening in & s transitions, we shall now propose a partial solution for 
these RGE's, partial in the sense that we will not provide a general solution for operators of dimension 6 but 
only recover the anomalous-dimension matrix for those operators that are usually considered. On the contrary, 
for dimension 7 operators, we will present the full result. This section shall also provide us with the opportunity 
to discuss the limits of our approach and how our extended analysis may be combined with the far-more advanced 
one of dimcnsion-6 vector operators. 



3.1 Dimension 6 operators 

First of all, we wish to check whether we can recover the usual leading-order anomalous-dimension matrix for 
dimension 6 vector operators. We thus consider the basis proposed in Eq.(2.15) of 12 . These can be viewed as 
specific linear combinations of ours (including also a rescaling) , so that we should be able to read the corresponding 
anomalous-dimension matrix 7^ from our results. For the time being, we forget about the dimension 7 entries. 
For the four-quark operators Oi-e of [l2], we obtain: 
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(26) 



which coincides with the first 6x6 subblock of Eq.(3.4) of this same reference. 

Similarly, we consider the so-called electroweak four-quark operators Q7-10 of Eq.(VII.2) of |3 and recover the 
matrix elements presented in Table XIV of this same reference. 

Let us turn to the magnetic and chromo-magnetic operators. In our approach, they satisfy the RGE, with 

{CG) = iCE,CQy 



\L.B.. 



diCc) _ 
din as 


^2(3) 1 

/3o 1 


(Cg) 


" 5 


0" 

3 
4. 
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. 3 
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- (C«)X«9 



71 



16 

-31 



(27) 



where we have neglected ^ ^ 1 for q b and replaced Vp^^^^ explicitly. The contributions from the 4-quark 
operators that we obtain here have no equivalent in the usual approach: the corresponding operators are simply 
not considered in the classical case, so that it is pointless for us to keep track of those hercp^ Therefore: 



^-(A^) = (^) 



5C2(3) 



3C2 (3) 
4^0 



Cq(mo) + (...) 



(28) 



Cq(mo) + (...) 



^*They are s ymboli cally replaced by (...) in Eq.| |28| l. Note that A was defined in the previous subsection and appears through 
^Fier^- see Eqs.(2l|25l. 
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To recover the usual scaling of the operators O7 
to the additional factor ( — 



[121, one has to multiply those coefficients by — , which leads 



^s(m) 
13(^0) 



11C2(3) 
4^0 



This is, again, consistent with Eq.(3.4) of 



12 



for the subblock 7 - 



Yet, one notices at once an apparent discrepancy with Eq.(3.4) of 12 : operators Oi_6 do not lead to divergences in 
the Ot^s directions. This is not surprising though, in the sense that the corresponding off-diagonal elements appear 
at the two-loop QCD level while we considered only the one-loop divergences. One may wonder why such two-loop 
effects in p2] are competitive with the diagonal one-loop elements of the O-^-g subblock and, if so, whether it 
endangers the validity of our one-loop analysis. The explanation is to be found in the relative order of the SM 
matching for the four-quark and magnetic operators (considering our normalization in Eq.(5[6l): if both were of 
the same order, then the two-loop mixing would only generate a subleading effect, of relative importance 0{as), 
which, for consistency, could be neglected. However, if one refers to the magnitude of the matching elements that 
we sketch at the end of section [2?Tj which is relevant for e.g. the SM, one observes that the matching coefficients 
of the magnetic operators are already of order as (translating the fact that it arises at the loop level even though 
the as factor itself is only an artefact of the relative normalization of the operators) , so that the two- loop mixing 
effect, related to the 0{a^) four-quark matching, becomes competitive: in other words, the naive hierarchy has 
been destabilized by the SM matching. More pragmatically, this situation in the SM is related to the fact that the 
closure of the charm-quark loop does not entail any as suppression. In our naive approach, we would naturally 
miss such an effect. Nevertheless, one also observes from the estimates at the end of sectio n|2.1| that this unbalance 
among matching conditions does only concern the dimension 6 four-quark operators Oi_ at least in a SM-like 
model: the corresponding matrix elements are therefore already known and may be included straightforwardly 
within our analysis. This potentially misleading aspect of our naive approach should, however, be kept in mind in 
order to combine our findings consistently with the usual dimension 6 analysis, directed at realistic high-energy 
models. 

Let us finally consider the semi-leptonic operators (C') — {Cg, Cy,Cip). They satisfy the RGB: 



rf(g') ^ ^2(3) 
din as /3o 



leading to: 
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(29) 



(30) 



.) stands for the inhomogeneous terms due to four-quark operators which we discuss later. If one factors out 
(for 5') or ag^mt (for F',T') in the definition of the operators, the scaling of the coefficients is modified 

_23C2jJ3) 11C2(3) 
^_s{p)_\ ""0 g^jjjj { as{p}_\ ^^0 



by factors 



/ as{p)_\ 



/ as 

^sit^o) J "''^^'^ V as (m o) J 

scalar and vector coefficients: see e.g. |30] 

O9 = (bj^PLs){h^l) and Oio = {brPLs){h^75l)- 

dimension matrix, translating the mixing with operators Oi-g, are given by 7^ = ~i' ~i' IT' ~i' "gJ ' 

7^ = 0, i = 1, . . . , 6. Up to a factor 2 accounting for the relative normalization of the operators, these matrix 
elements coincide with Eq. (VIII. 11) of [s]. 

We have thus verified that, except for the mixing of four-quark vector operators with the magnetic and chromo- 
magnetic operators, which is a two-loop effect, our calculation covered consistently all the leading-order elements 
of the anomalous-dimension matrix for those operators that had been considered in the literature. 
Before turning to the solution for dimension 7 running, it is worthwhile discussing how dimension 7 corrections to 
these dimension 6 RGE's should be implemented. First note that we are talking here about corrections of order 
at most, so that it would make little sense to include them if the dimension 6 coefficients were considered at 
leading order only: dimension 7 effects will make sense numerically only if dimension 6 contributions are known 
up to 0{as) or even 0{a'g) (depending on the matching conditions). A second remark comments on the absence 
of dimension 6 effects in the RGE's of dimension 7 operators (in other words, the anomalous dimension matrix 
is block-triangular): this naturally proceeds from the fact that the dimension 6 basis is stable and does not 



(respectively), leading to the commonly known non- running of 

Let us now consider the inhomogeneous terms in the classical case of 
Then we find that the off-diagonal elements of the anomalous- 

3, are given by 7^ = iJ, [- ^ ' ' ' 



Actually O2 only, at the level of the matching conditions, but, O2 closing on the whole Oi-e 
whole subset needs to be included. 



basis under renormalization. the 
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require dimension 7 counterterms to cancel its divergences (this can be derived from simple power-counting) . The 
dimension 7 RGB's can therefore be considered separately while their mixing-effect is injected directly as a small 
inhomogeneous term in the dimension 6 RGB's. 

3.2 Dimension 7 RGE's 

We present the solution for the RGB's of dimension 7 operators. We stress again that those are independent from 
the dimension 6 running, allowing us to solve them separately. The block-triangular shape appearing in 7^'" the 
restriction of the anomalous dimension matrix to dimension 7 operators, invites for a splitting into subblock RGB's 
with inhomogeneous terms: this observation makes an explicit solution tractable analytically, which we perform 
for the sake of completeness. Note that the methodology, consisting in splitting the anomalous-dimension matrix 
according to its block triangular shape and resulting in inhomogeneous linear equations, is a recurring feature in 
BBT's and has been used, e.g. in [3l] (although in a different context; see section 3.4.2 of this reference). 

3.2.1 (fes)-Gauge operators 



Gluonic operators] (Cg) = [Cqd ,Cqd ,Cqs ,C qs tCqt)^'^: 
The RGB is a simple homogeneous Ist-order differential equation: 
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(31) 



with Aqq = OqVqOq ^ where Vq is a diagonal matrix: 

Vq = diag [-^t, 1(117 - ^/3073)l, ^(277 - \/3T93)l, ^(277 + V'STM)!, i(117 -I- V3073)l] 
The details of this diagonalization can be found in Appendix [C] 

Hybrid operators; {C-u) = (C^s, C^^s, C^r)-^'^: 
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The exponentiation of A-un can be achieved without difficulty through its diagonalization: Ann — OuDuO-j^ , 
with V-H = diag [-^l, |l, ^l] (refer to Appendix [c]) . 
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Interestingly, certain of those operators have a negative anomalous dimension, meaning that they will grow at 
low energy, and thus become important. Note however that this negative dimension may simply be an artefact 
of artificial factors a^^ in the normalization of the operators: the matching conditions in a definite model are 
therefore likely to correct this feature by introducing suppression factors. Otherwise, such growing directions in 
the dimension 7 anomalous-dimension matrix could lead to significant perturbation of the dimension 6 results. 
While the computation of the anomalous dimension is particularly simple, e.g. in the case of the photonic operators 
(leaving little room for false moves), the prospect of dimension 7 contributions enhanced at low energy should be 
considered however with caution. 

3.2.2 Dimension 7 (bs){qq) operators 
fr<i \ — (ni ri ri ni \l,rl,r. 
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Aqq can be explicitly diagonalized: Aqq — C'qdiag(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, |l?24)Cg ^ (refer to Appendix 
ICl). 



3.2.3 Dimension 7 (bs){ll) operators 
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Note that the positive eigenvalues (negative anomalous dimensions) change sign once the artificial factor a^^ is 
factored away from the operators. 

This concludes the solution of all dimension 7 RGB's. We will now consider a simple example engaging these 
dimension 7 effects. 



4 A simple application 
4.1 Setup 

Let us assume that leading new-physics effects would arise in the possible (though unlikely) form of the following 



dimension 6 SM operator 27 



Ml 



(39) 



where g' and g are the f/(l)y and SU{2)l coupling constants, Bf^^, and W^^t"- the corresponding field-strength 
tensors (r° correspond to the SU{2)l generators), while H is the Higgs doublet (with negative hypercharge) . Kgi, 
is a complex coefficient encoding new-physics effects. At low-energy, it results in a modified Z — b ~ s coupling: 



(40) 



where v = (2a/2Gf)^^^^ is the electroweak vacuum expectation value. Such a term could generate flavour effects 
at the Z-pole, resulting in a perturbation of the precision measurements [Z — >■ bs, sb being mistaken experimentally 
for a Z — >• bb decay). 

Considering the modified Z-coupling, we compute the associated matching conditions, at tree-level, for the oper- 
ators of the b s transition: 
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(41) 

where we here use the following normalization of the operators (/ stands for any of the low-energy quarks and 
leptons): 



(42) 



Note, at this point, that the classical framework of dimension 6 vector- type operators would not allow for a 
consistent description of these new-physics effects: albeit the (y-'')^'^ operators receive a contribution, the 
associated effect is of the same order as that of the dimension 7 (H^)'"'^ operators. One should thus rely on 
our extended analysis. 



4.2 BR{Bs 

Now let us consider the decay Bg — > l^l^. This rate is one of the traditional search channels for new-physics. 
Evidence for its observation (in the case / = /i) has been reported a few months ago at LHCb [l], evidencing a 
good agreement with the SM, hence constraining new-physics effects more tightly. For a recent review, refer to j2]. 
The well-known SM-matching provides [sj (the function Y is defined in this reference): 



(Mz) 



ITS 



w 



-6.2 • 10" 



(43) 



We (safely) neglect dimension 7 effects associated with the SM and assume that only new physics associated with 
Ksb may lead to a significant deviation. 
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Then, we consider the foUowing matrix elements (in terms of the Bg-meson decay constant fs, and its four- 
momentum P,,): 



Ml 



{0\h^s\Bs) = -iJb- 
(0| 67^755 =z/bP^ 
(0| h{t - ^)^75S \Bs) - ifBinib - m,)P^ 



\bs\B,)=Q 
I h^s \B,) = 



(44) 



32 



and derived the resuh for the matrix elements (0| 6i(i5— S)^1/75S \ Bs 



Here we have followed the conventions of 
One then derives (where we consider only the operators and i?', relevant in this particular example): 



BR{B, ^ m-) ^ '^^'"^ " ^ ~ ^^''^ ^ ^^'"^ ^ 



LL 



Note that for the combinations C, 



L.RL 

v 



(45) 

Cyf - and C^f^ - C^f, the mixing to the four-quark operators 



RR 



cancels in the RGE (Eq.( 16 1, due to the 1 + S in TW'^') so that we obtain simple scalings for both these quantities: 

{^)' = 1 for C^' - and (^)"'^ = {^Y^' ^^.^ " C'i.r ■ Note however the 

off-diagonal term in M-^^ mixing and operators. The RGE for the vector coefficients hence reads: 



a m as 



L,R.R^ 



2C2(3 ) 
/3o 



LMR-i 



l<-yi -<-v; \{^^b) = [Cy\ -Cy, \[Mz) T, — mb{Mz)[C^, ~C^\ \(Mz) 



In 



= [C^f^ - Ci;f'\Mz) - mbiMz)[C^f'' - C^f'KMz) 



^L.RR 



L,RL 



yL,RR^ 



»S(Mz); 2C2(3) ^ 

e dec 



2C2(3) 

/ Qs(Mfc) \ 

Us(Mz) 



- 1 



(46) 



3C2(3) 



It is also more convenient to work with a low-energy b mass: mb — mb{Mz) ^ ^sImI) ) ■ ^o^e that this 
replacement does not concern the explicit factor rab multiplying the iJ' coefficients in Eq.(45). Replacing the 



coefficients at the Mz scale by the matching conditions, we then observe a complete cancellation of the new- 
physics contributions: 

flMBm ' ~ 



BR{Bs r r) 



327rrB 



(47) 



The scenario under investigation hence receives no constraint from BR{Bs — > l^l^). Although this is what the 
naive tree-level calculation would have predicted, it is a non-trivial result to observe that this feature is preserved 
by the resummation of the leading logarithms through the RGB's: inconsistently neglecting the H contributions 
or the V — H mixing, for instance, would have generated limits of the order \Ksb\ < 10^"^. Note however that 



the cancellation of the new-physics effects is tightly related to the form of the matching conditions of Eq. ( 41 1 and 



that a perturbation of Eq.(41), e.g. through the implementation of additional dimension 6 operators, would lead 
to relevant limits again. 



4.3 BR{B Kui>) 

Many other b s transitions are of course relevant to constrain Kgb- Let us consider another simple example: 
the decay B Kvi? is bounded by the limit BR{B- K' + Inv.) < 13 • 10"^ (90% C.L.) from BABAR [33]. 



The SM provides BR{B- R-vv) ~ (4.4j:i:?t°;^t°;^) • 10"^: refer to 
coefficient for the SM can be read in |i3 (the function X is defined in t 



34 



for a summary. The relevant matching 



lis reference): 



C^HMz) 



%/2G 



us 



pa ^ , m 
VtlVtbXi 



w 



1 ■ IQ- 



(48) 
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Concerning the decay constants, we again follow the conventions of 32 and derive the matrix element correspond- 
ing to the _ff-type operators: 

A/f2 _ 

< K{PK)\s^^b\B{PB) >= [Pb + Pk]M^) + HPb - Pk]Ms) - f+{s)) s = {Pb - Pk? 



< K{Pk)\s-<j,,,{Pb - PKYb\B{PB) >= i<[Pb + Pk 



< K{Pk)\si(J^ - ^)^b\BiPB) >~ [Pb + Pa'Ip | mbf+{s 

Ml- 



Mb + Mk 



[PB-PKUMB~MK)>fTis) 



Mb + M 



K 



(49) 



[Pb - Pk], 



rrib- 



K 



{his) - f+{s)) + {Mb - Mk)It{s) 



For the form factors, we follow the discussion in 34 closely, employing the paramctrization of [35] and the 
approximate relation /t(s)//+(s) ~ ^^^^g^- 

The estimation of the branching ratio is now straightforward {Ny = 3 is the number of neutrino flavours) : 



BR{B Kvv) = 



bl2^ni%TB Jo 



2 a2i3/2 



ds [s^ - 2s(M| + Ml) + (Ml - M|.)2] 



■sM-s) 



iC^i^ + Cfji^) 



Mf 



K 



(Aifc) (50) 
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where we will use Mg = 5.27925(17) GeV, tb = (1.641 ± 0.008) • lO^^^ g, A% = 0.493677(16) GeV 
The running of the coefficients is again very simple as the neutrality of the neutrinos ensures no mixing with the 
four-quark operators. Similarly to Eq.(46l, we have: 



[Ci;i^ + C^,/^](m6) = [Ci^i- + C^^Mz) - m,{Mz)[C'^t + C^i^KMz) 



as{f^b) 
as{Mz) 



2C2(3) 



[Cji!^ + cSi^]{^,, 



C2(3) 

\as{Mz) 



[C^^^+CmMz 



(51) 



Going to our explicit matching conditions of Eq.(41 ), we observe that new-physics effects here persist in the part 
multiplying /r(s). Numerically, we find Ksb (which we assume to be real) in the range [—4, 1.5] • 10"'^. The 

naive estimate Kgb ^ 10~^ ^ returns Ajvp ^ TeV; this scale can be lowered further if one assumes that 

the new-physics operators should be loop-suppressed, follow Minimal Flavour Violation, etc. The bound on new 
physics appearing in this fashion is therefore relatively loose. Note however that several other observables in the 
5 — s sector should be considered before drawing any conclusion concerning the viability of Ksb ~ 10"'^. Since 
our focus in this section was merely to consider a concrete, yet simple, case where the inclusion of the dimension 
7 RGB's was relevant, we shall not pursue this analysis further. 



Let us briefly summarize our achievements, we have established the most general basis of (on-shell) operators, for 
the EFT describing 6 — > s transitions, up to mass-dimension 7. After computing the associated ultraviolet QCD 
divergences at leading order, we have derived the corresponding RGE's. Comparison with the existing studies 
concerning dimension 6 operators proved satisfactory. We have also solved all the RGE's describing the evolution 
of pure dimension 7 operators: interestingly, we observed that some directions exhibited negative anomalous 
dimensions, which could lead to an enhancement effect at low energy. Note however that this property depends 
on the normalization of the operators so that the matching conditions in a definite high-energy model are likely 
to regulate it. We finally used this analysis to constrain, using the measurement of the decay 5° — ^ A'^Ai^ at 
LHCb and the limit set by BABAR on S — > Kvv, a very naive extension of the SM resulting from the addition 
of a dimension 6 SM-operator. More generally, let us recall that the inclusion of dimension 7 effects has little 
relevance in e.g. the SM, where such operators receive extra-suppression due to the very-constrained pattern of 
flavour- violation. Beyond the SM, requiring New-Physics to project only on operators of higher-mass dimension 
would be an elegant way to circumvent the strong limits on non-standard flavour violation: however, such models 
would have to be designed on purpose. 
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A Fierz identities 

We derive/recover the following Fierz identitie^^ 
• Scalar identities: 



ii^'Pn^m'PRX) = U^'PRxMPni^) - 3^(V^V^-Pflx)(xV^,PflV) 
{^'Pl^WPrx) = Wl''PRX){x'lt.PL^) 
WPr^){x'Plx) = h(^'l''PLX){x'l^PRi') 



(52) 



Vector identities: 



'{i^'l^PLi^){x'lt.PLX) - -(V5V^lx)(x'7m^lV') 
{i^'l>'PR'iP){x'li.PRX) = '{i''l''PRX){x'l^.PR^) 
ii^'rPL^Wlt.PRX) = 2{i:'Pnx){x'PLi^) 

X^'i^Pr^){x!i,.Plx) = 2{^'PLX)ix'PRi^) 



(53) 



Tensor identities: 



(V;V''-Pz,V)(xV^.Plx) ^ 
{^'a''''PBMx''J^.PRX) 



-2^{fPLx){x'PL^) - Wa>'"PLX){x''^t..PL^) 
-24(V7Pflx)(x'^i?V') - W^'"'PRX){)^'^t..PR^) 



(54) 



Hybrid identities: 



{i^'PLiD'^r^)(^-f^PLx) = U^'PLxKx'PR^Pi^) - -i^'a^^^PLXm'l.PLiD^^) 
{^'PRiDt^^){x'i^,PRX) = \{^'Prx){)^PliP^) - -^{^'''^''''PRxWluPRiD^i^) 
{ij'PLtD'^^){x'j^PRX) = ii''rPRXWPLtD^^) - ^^i^'rPRX)ix'-f^.PR^Pi^) 
[{i>'PRiD^^){x'j^PLX) = {^''rPLxKx'PR^D^.^) - Wl^'PLXm'lt.PLip^) 



(55) 



One can straightforwardly generalize these results to the case where the covariant derivative acts on the first 
spinor, instead of the second. Moreover, with a little bit of algebra: 

{^'(y^''PL.,RX){x:i.PL,R^D,,i^) = - l'5)^PL,RXKx'l^PL.,R.i^) - - iS)M75PL,/?x](xV75PL,K^) 

- [^'(^^7^PL,i^ - i''Pr.li^)x]{x1^Pl,r^) " W{^%Vlm + i''ToPr.li~^)x]{x lt.l^PL,R^) 



l{^y'J^''PL,RX)[x'{%^.PL,R 



- ^^..PR,LP)i^] - ::^^p[{^^ya,,,PL,Rx){x'l''Yh''PL.R^] (56) 



Consequently, for q — b,s, 0\ ^ ^ ^ = 5, V, T, H^^fl l^r can be expressed in terms of 0\ r l r = S,V, T, H, 
ji jf sls, fi ~ "^ierz^i- R l r (*^^^ ^^S^ originates from the necessity of anticommuting fermion fields), 



-'^Note that, here, we consider the identities only from the point of view of the Clifford algebra, without including the — sign resulting 
from the anticommutation of fermion fields. 
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and reciprocally. Chiralities are ordered as LL, LR, RR, RL. 
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(58) 



B Anomalous Dimension Matrix 

Considering that RGB's in EFT's already represent a significantly equiped industry, we have decided to collect 
here our results for the anomalous dimension matrix 7c so that one may easily implement this matrix without 
having to refer to our discussion in section [2]3]2] and sequels, oriented at a more progressive, if exhaustive, solution. 
Note however that, due to the size of the matrix, we are bound to present it block after block. Remember that 
our operators are ordered as in Eq.( 4|5|6 1 and that the chiralities follow the ordering LL, LR, RR, RL, or L, R 
(depending on the number of chirality indices the operator carries). The definition of the blocks in chirality space 
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is reminded at the very end of this section. RecalUng Eq.(15): 










" 


ic = 


ql 

Tc 




7c 




gi 
lie 




QQ 

Tci 



(59) 



• Subblock 7^ 



We remind the reader that operators with different lepton flavours do not mix (at least as long as massless neutrinos 
are considered) so that the subblock 7" is diagonal in lepton-flavour space: 



7c 



7c 

















r' 

















Ic 


































7c 
























7c^^ 



(60) 



From Eq.(16) and sequels, we read, for any lepton flavour I (but note that the basis of operators is shortened in 
the neutrino case, due to the absence of low-energy right-handed neutrinos): 



2asC2{3) 
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fl 
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2(1 + ^^) 

\ nib J 
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(61) 



Subblock 7j, 



nl 

Tc 



7c 

Ic 
.Ic 



~ Ull 

7c 

~ da 

7c 

- SLL 

7c 

~ CLL 

Ic 

~ ba 

ic 



Tc 
lis 

Ic 
Ic 



Ic 
Ic' 

Ip 

Ic' 
Ic 



Ic 

Ip 

Ic' 
Ic' 
Ic' 



Tc 

Ic^ 

Ic^ 

Ic^ 

Ic^ 



(62) 



From Eq.( 16 ) and sequels, we read, for any lepton flavour / and q = u,d, c: 
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2qsC2(3)3 
47r 2 



QlQq 



"0 













0' 
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(63) 



For q = b,s, we have: 



~ql 20562(3)3 

7c = 4^ 2^iQg { 
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(64) 
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Subblock 
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7c = [ic Ic" Ic 



Tc 



From Eq.(16) and sequels, we read, for any lepton flavour /: 
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7c 
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Tc = 
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7c 

~du 
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~su 

7c 
Ic 
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~dd 

7c 

~sd 

7c 
Ic 
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Tc 
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Tc Tc 



Tc 
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For two quark flavours q and q' , we read from Eq.(20l and sequels, 
* \i q and q' ^ b, s: 



Ic 
Tc' 
Tc' 



7c 



47r 



{(5'?'«diag(51, 51, 51, 21, 21, 51, 51, 51, 21, 21) - ^^'^TW^'f^ - A^lt. } 



if q and q' = b ov s: 
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5''«[1,0]A^L'« 
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if g 7^ 6, s and q' ^ b or s: 



if q — b or s and g' 7^ 6, s: 
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• Subblock j^? 



- qd 

Tc 



- as 

Tc 



~ gc 

Tc 



For a quark flavour q = u,d, s, we read from Eq.(20) and sequels: 
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For q ^ b, s, we have: 
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Subblock 7^^ 
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For a quark flavour q = u,d, s, we read from Eq.(23) and sequels: 



For q = b,s, we have: 
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We finally collect the definition of the various relevant subblocks in chirality space (with exception of the identity 
1 and the null matrix 0): 
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C Diagonalization of the anomalous-dimension matrix for dimension 
7 operators 

C.l Gluonic operators 



We consider the matrix Aqq defined in Eq.(31 1 
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Its diagonalization happens to be accidentally tractable 
^1 i(117-N/3073)l,;i(277 
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so that the exponentiation is fully determined. 
C.2 Hybrid (6s)-Gauge operators 



We consider the matrix A-hv. defined in Eq.(32) 
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Its diagonalization is straigtforward: 
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allowing for a simple exponentiation: 
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as(Aio) 



5C2(3) 



coshi?(/i)l 
—2 sinhi?(/i)cr3 





I sinhi?(/i)cr3 
coshz9(/i)lL 







6102(3) 
24/^0 



^(,).^^lni^ 



C.3 {bs){qq) sector 



We consider the matrix Aqq defined in Eq. ( 36 1 : 



An 



FierzJ 



Let us introduce the following blocks: 



+ 1 + S 



^Fierz 



q—b,s 



C',diag(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2?24)C<7^ 



























-1 





1 





















1 










n = 


1 


1 





1 
















~ 


















1 


-1 



























1 











1 


1 












1 








-1 





1 












-1 














-4 








-1 














-4 












1 








-1 





1 











-1 


1 














1 












-1 














-4 








-1 














-4 




1 





-1 


1 














1 














1 













71 


-1 








-1 

















-1 


-1 





































-1 























1 












-1 


-1 

















-1 








-1 






























-1 












31 



1 


























1 




















1 

3 





n 


























1 




























































































































1 


















































1 












































1 

3 





n 


















































1 


















































1 


















































1 












































1 

3 





n 


















































1 


































































































w 







































































































1 























































1 























































1 



































2?24 





^24 













1 



























































1 























































1 




















































1 



















2?24 





P24 





























1 

























































1 



















































1 



















































1 


-V2i 


2?24 





^24 










































1 



















































l-¥2?24 

















































1 


I?24 







































-2?24 





-fl?24 





1 



















































1 



(87) 



Oq turns out to diagonalize Aqq satisfactorily. 
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